Long-range crystalline nature of the skyrmion lattice in MnSi 
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We report small angle neutron scattering of the skyrmion lattice in MnSi using an experimental 
set-up that minimizes the effects of demagnetizing fields and double scattering. Under these condi- 
tions the skyrmion lattice displays resolution-limited Gaussian rocking scans that correspond to a 
magnetic correlation length in excess of several hundred /xm. This is consistent with exceptionally 
well-defined long-range order. We further establish the existence of higher-order scattering, discrim- 
inating parasitic double-scattering with Renninger scans. The field and temperature dependence of 
the higher-order scattering arises from an interference effect. It is characteristic for the long-range 
crystalline nature of the skyrmion lattice as shown by simple mean field calculations. 

PACS numbers: 75.25-j, 75.50.-y, 75.10-b 



Single-crystal small angle neutron scattering (SANS) 
in the A-phase of the itinerant helimagnet MnSi recently 
established a highly unusual six-fold symmetry of the 
scattering pattern perpendicular to the magnetic field 
B [1 . Regardless of the orientation of the crystal lat- 
tice with respect to the magnetic field the same six-fold 
diffraction pattern was seen, characteristic of a magnetic 
structure which is almost completely decoupled from the 
underlying atomic lattice. A theoretical calculation in 
turn identified the A-phase in MnSi as a skyrmion lat- 
tice stabilized by thermal fluctuations, i.e., a new form of 
magnetic order composed of topologically stable knots in 
the spin structure. The associated non-trivial topology 
was confirmed by means of the topological Hall signal [2] . 
Both the theoretical analysis [1] and small angle neutron 
scattering in Mni-a^Fe^^Si, Mni_a;Coa^Si and the strongly 
doped semiconductor Fei-a^Coa^Si further suggested that 
the skyrmion lattice is a general phenomenon [3], with 
spin torque effects at ultra-low current densities as the 
most spectacular result [4]. 

These studies inspired Lorentz force microscopy pro- 
viding direct evidence of skyrmion lattices and individ- 
ual skyrmions in Fei_a;Coa;Si [5 , FeGe [6 and MnSi [7 . 
In the thin samples and for the perpendicular magnetic 
fields studied, the skyrmion lattice is thereby more stable 
^ than in the bulk, but with increasing sample thickness 
the magnetic phase diagram was found to approach that 
of bulk samples. This suggests that the same skyrmion 
lattice is realized in bulk samples and thin films. 

Taken together a pressing question regarding the exis- 
tence of skyrmion lattices in bulk materials concerns di- 
rect microscopic evidence of their long-range crystalline 
nature, i.e., the precise spatial variation of the magne- 
tization on very long scales. The quality of the long- 
range order strongly influences both the pinning forces 



[H [8] and the rotational spin torque forces observed in 
the presence of currents [4] . 

In principle, quantitative information on the magnetic 
structure can be obtained from a reconstruction from 
higher-order peaks in neutron scattering. However, pre- 
vious neutron studies [1 were subject to strong dou- 
ble scattering, i.e., neutrons scattering twice from the 
magnetic structure, which are difficult to discern from 
higher-order peaks. Moreover, an exponential (rather 
than Gaussian) intensity variation of the rocking scans 
and streaks of intensity emanating radially from the first 
order peaks [1 seemed to hint at unusual aspects of the 
morphology of the spin structure. 

In this Letter we report a high-resolution small angle 
neutron scattering study of pure MnSi to resolve these 
issues. As the skyrmion lattice is extremely weakly cou- 
pled to the atomic lattice and follows closely the applied 
magnetic field, it is essential to guarantee a homogeneous 
magnetic field inside the sample by reducing all effects of 
demagnetization fields. We used therefore a thin sample 
with a thickness of '^ 1 mm which was illuminated in a 
small central section only. This also allowed us to reduce 
the amount of double scattering. Taking these precau- 
tions we find sharp resolution-limited Gaussian rocking 
scans, while all peculiar features such as streaks of inten- 
sity vanished. Our results highlight, that extreme care 
has to be exercised before claiming an unusual morphol- 
ogy of the skyrmion lattice in any B20 compound. 

As our main result we unambiguously establish higher- 
order scattering, by discriminating parasitic doubly scat- 
tering in Renninger scans. The magnetic field B and 
temperature T dependence of the higher-order scattering 
may thereby be quantitatively explained with an interfer- 
ence effect, providing unambiguous microscopic evidence 
of the long-range crystalline nature of the skyrmion lat- 
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FIG. 1: Typical SANS data in the A-phase of MnSi. (a) In- 
tensity pattern in a thick cyhndrical sample at fioH = 0.16 T 
as reported in [l] . (b) Rocking scan of the data shown in panel 
(a) as reported in [l]. (c) Intensity pattern in thin platelet 
(sample A) at fioH — 0.2 T. (d) Rocking scan of the data 
shown in panel (c); note the narrow Gaussian dependence. 



tice in bulk samples of MnSi. 

For our studies two MnSi platelets denoted as sample 
A and B were cut from the same ingot used previously 
p^miOHTT]. Sample A was ~ 14 x 9 x 1.4 mm^, and sam- 
ple B was ^ 12 X 7 X 1 mm^, both with a crystalline (110) 
direction normal to the platelet. The crystalline mosaic 
spread was measured to be very small, ~ 0.15°. The 
specific heat and resistivity of test pieces from the same 
ingot are in excellent agreement with literature with a 
good residual resistivity ratio, RRR ~100. The samples 
were prealigned by x-ray Laue backscattering. For pre- 
cise alignment in the SANS studies at low T we use that 
the magnetic structure at zero field is described by helices 
in (111) directions. 

Our studies were carried out on the cold diffractome- 
ter MIRA at FRM II (Munich) and the SANS instru- 
ment V4 at BENSC (Berlin). Neutrons with a wave- 
length A = 9.7 A ± 5% and A = 4.5 A ± 10% were 
used at MIRA and V4, respectively. The instrumen- 
tal resolution at MIRA was A/^^^ = 4° in azimuthal 
direction, A/3q = 0.004 A in radial |q| direction and 
A/3kf = 0.35° perpendicular to |q| in the direction of 
kf. This compares with a resolution of /S.j3az = 4.9°, 
A/3q = 0.003 A-\ and A/3kf = 0.21° at V4. Tempera- 
tures are given with respect to the helimagnetic transi- 
tion temperature Tc determined by SANS. 

The importance of the sample thickness for SANS 
studies is evident from previous work in the helical 
state, where distinct double scattering was observed, e.g.. 
Fig. 2 (A) & (D) in Ref. p. Consistent with the entire 
literature on SANS in the zero field helical state (cf. 
Ref. [121 [E] and references therein), typical data in the 
helical state of sample A (not shown) display a Gaussian 
rocking dependence with a rocking width of rjrn = 1-6° 
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FIG. 2: Operating principle and typical data of Renninger 
scans, (a) Ewald-sphere depiction of Renninger scans, see text 
for details, (b) Typical scattering pattern obtained by a sum 
over a rocking scan around (j) after background subtraction 
recorded at high T. (c) Intensity as a function of rocking 
angle in a Renninger scan. The intensity was integrated 
over the areas indicated by box 1 and box 2 in panel (b). 



corresponding a correlation length ^ 10^ A. However, 
due to its reduced thickness double scattering is strongly 
reduced in sample A. 

The effects of the sample shape are most dramatic for 
for SANS studies of the A phase as illustrated in Fig.[l](a 
& b), which reproduces Fig. 2 (E) of Ref. [1 . These pre- 
vious results suggested a rocking width similar to the 
helical state, with an unusual non-Gaussian line shape 
(Fig.[l](b)). For another sample studied in Ref. [1 strong 
double scattering and even unexplained streaks of inten- 
sity emanating radially from the first order scattering 
were seen. In stark contrast, we do not find any such 
peculiar features for our thin platelets (see Fig.Il](c & 
d)). Moreover, the rocking dependence is now Gaussian 
with an extremely narrow width, tja = 0.45°, slightly 
larger than the resolution limit A/3k/ = 0.35°. Thus, 
the intrinsic magnetic correlation length of the skyrmion 
lattice exceeds 100 /im and is therefore more than a fac- 
tor of 100 larger than for the helical state. In contrast, 
previous studies reflected the variation of the internal 
magnetic field directions due to demagnetizing effects. 

To distinguish double scattering from higher-order 
scattering we used so-called Renninger scans depicted in 
Fig. [2|a) [14]. The sample is thereby first rotated to- 
gether with the magnetic field around the vertical axis 
through an angle x until the sum of two scattering vec- 
tors qi + q2 touches the Ewald sphere thus satisfying 
the scattering condition. This is followed by the actual 
Renninger scan, which is a rocking scan with respect to 
qi + q2 through the angle ^, while recording the inten- 




FIG. 3: Typical Renninger scans at various T and B in the A- 
phase of MnSi. Panels (a) through (c) were recorded at T = 
Tc - 1.5 K and magnetic fields /j^qH = 180 mT, 200 mT and 
220 mT, respectively. Panels (d) through (f) show Renninger 
scans for T = Tc — 0.5 K and magnetic fields of /hqH = 180 mT, 
200 mT and 220 mT. True higher-order scattering is shaded 
in red and remains constant for large 0. 



sity at qi + q2. This way double scattering is "rocked 
out" of the scattering condition, while higher-order scat- 
tering continues to satisfy the scattering condition for 
all 0. For the Renninger scans sample B was mounted 
with its crystalline (110) direction parallel to the incom- 
ing neutron beam. A crystalline (100) direction was ori- 
ented approximately horizontally. The background was 
determined for T well above Tc for each rocking angle and 
subsequently subtracted. The intensity at qi+q2 as indi- 
cated by box 1 in Fig.^panel (b) was then compared with 
the intensity in a box of equal size at a position slightly 
to the side of qi + q2, labelled box 2. Typical varia- 
tions of the intensities in box 1 and box 2 with the angle 
(j) are shown in Fig. [21 panel (c) for T = Tc — 0.5 K and 
jhqH = 200 mT. The intensity observed at qi +q2 clearly 
displays two contributions: (a) two Gaussian peaks due 
to double scattering when either qi or q2 intersect the 
Ewald sphere, and (b) a constant intensity arising due to 
true higher-order reflections (red shading). 

Typical Renninger scans for T = Tc- 0.5 K, T^ - 1.0 K, 
and Tc — 1.5 K under magnetic fields of jhqH = 180 mT, 
200 mT and 220 mT are shown in Fig. |3] With increas- 
ing T and B the higher-order scattering increases. The T 
and B dependence of the scattering intensities are sum- 
marized in Fig.|4j Panel (a) shows the field dependence 
of the first order peaks at qi and T = Tc — 0.5 K. The 
peak height is shown, which was confirmed to be pro- 
portional to the integrated intensity. Gompared with the 
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FIG. 4: (a) Intensity of first order peaks in the A-phase as a 
function of B for T = Tc — 0.5 K. (b) Ratio of the intensity at 
Qi +q2 and first order diffraction at various temperatures, (c) 
Galculated staggered field at |bg^+q2 1^, Eq. (m), which induces 
the scattering at qi + q2 for different values of the phase a. 
(d) Galculated ratio of the higher-order diffraction at qi + q2 



and first order diffraction, 



|Mqi+q2l'/ 



iMa 



as a function 



of the magnetic field in mean field theory for t — —.9,-1 
and —1.1 in the units of Eq. ^\. Dashed (dot-dashed) line: 
|M2qjV|Mq,P (|M2q,+q,|V|Mqj2) fori = -1. 



field dependence for B || (111) reported in Ref. [15 for 
Tc — 0.2 K, we find a more gradual variation of the inten- 
sity. Panel (b) shows the ratio of the higher order peak 
intensity at qi + q2 with respect to the first order inten- 
sity at qi for T = Tc- 1.5 K, Tc - 1.0 K and Tc - 0.5 K. 
Above jiqH ^ 180 mT this ratio increases as a function 
of T and B. Unfortunately it was not possible to track 
the higher-order scattering below fi^H — 180 mT, since 
the first order intensity was too weak. 

In summary our main experimental results are: (i) A 
strong magnetic field dependence of the second order in- 
tensity, apparently vanishing for a certain field B ^ Bi^t 
inside the A-phase. (ii) An increase of the second order 
intensity with increasing T. This may appear counter- 
intuitive, since the non-linear effects leading to higher 
order peaks should be less pronounced when all ampli- 
tudes decrease with increasing T. And finally, (iii) a tiny 
weight of the higher-order peaks of the order of 10~^. 

For a theoretical analysis, we first discuss what can be 
inferred from the experiments without a detailed calcu- 
lation and compare in a second step our experimental 
results to mean field theory. Almost all of the scattering 
intensity arises from six resolution-limited main scatter- 
ing peaks, and therefore the magnetic structure is well 
approximated by a superposition of three helices and the 
uniform magnetization. 

A single pair of peaks at iqi describes a spin helix 
of a given chirality [16] determined by the chirality of 
the atomic structure. A representative of such a he- 



lix is -^(0,singx, — cosgx) or, in Fourier space, Mq^ = 
^(0, — z,— 1) for qi = (g,0,0), where <l>^ is the weight 
of the peak. Our observation of higher-order scatter- 
ing proves a crystahine character of the magnetic state, 
which may therefore be described by a hnear superposi- 
tion of three such hehces rotated by n27r/3 (n = 0, 1, 2) 
around the z axis. The phase relationship between the 
hehces thereby distinguishes the skyrmion lattice from 
other forms of magnetic order. Two of the three phases 
just describe translations of the magnetic structure in the 
xy plane. However, the third phase a, defined by mul- 
tiplying Mq^ by e*", while keeping the other two helices 
unmodified, strongly affects the magnetic structure (in 
our convention a = describes the skyrmion lattice [1 
where in the center of the skyrmion all three helices point 
antiparallel to the external magnetic field). 

Unfortunately, no information on the relative phase of 
the three helices can be inferred from first order scat- 
tering, since the signal is only sensitive to |Mq. p. In 
contrast, the higher-order peaks are very sensitive to the 
phase OL because they are subject to interference effects. 
When considering the Fourier transformation of the lead- 
ing magnetic interaction term M(r)^, collecting all terms 
linear in Mq^+q^, jM'^d^r = Mq^+q^bq^+q^ + ..., 
one obtains an oscillating effective field bq^+q^ arising 
from Mq^ g and the uniform magnetization Mq. From 
qi + q2 = 2qi + q3 = 2q2 + qe (see Fig. ^) one obtains 
that the interference of several processes determines the 
strength of bq^+q2. Adding all of these terms (examples 
are Mo(M_q^M_q2) or M_q^(M_ 



.q^Mqg)) we obtain 
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A plot of this function for various values of a is shown 
in Fig. [4](c). As negative Mq are unphysical, one finds 
only for values of a close to a = that bq^+q^ be- 
comes very small for a certain magnetic field where 
M^l^ = a/2 24/37 ^ 0.92 (to be compared to 0.94 and 
0.96 obtained from the mean field theory given below for 
t = — 1 and t = — 5, respectively). The observed strong 
suppression of Mq^+q2 is therefore the signature of a very 
small a characteristic for the skyrmion lattice [l]. 

A complete quantitative theory of the skyrmion phase 
requires account of the effects of thermal fiuctuations [1] . 
Yet, a semi-quantitative explanation of our experimental 
observations may already be obtained on the level of a 
simple mean-field approximation. After a rescaling of co- 
ordinates, magnetization, magnetic field and free energy, 
the Ginzburg Landau free energy density in the presence 
of the Dzyaloshinskii-Moriya interaction ~ M(VxM) is 
given by [H HH Hg 

F = (l + t)M2 + (VM)^ + 2M(VxM) + M'^-BM (2) 

where t measures the distance to the B = mean-field 



critical temperature. Other contributions can be ne- 
glected as they are higher order in spin-orbit coupling or 
give only small contributions close to the critical temper- 
ature where the skyrmion lattice is stable. Minimizing F 
with the Ansatz M(x) = En,m e'^'^'''^"''''^" M^qi+mq2 
for integer n and m provides the relative weight 
of |Mq^+q2p/|MqJ^ aud othcr higher-order peaks, cf 
Fig.[4]^d). An open issue for future studies is the gradual 
variation of the amplitude Mq^^ for small 5, cf 



Fig.ma, 
on [fl. 



where theory predicts a sharp first-order transition 

Our calculation reproduces the main experimental ob- 
servations, including the tiny weight of the higher-order 
peaks and the approximate vanishing of the signal for a 
certain magnetic field, B ^ ^int due to the above dis- 
cussed interference effect. Most importantly, it explains 
the unexpected rise of the signal for increasing T: the 
overall drop of the higher-order scattering for larger T is 
overcompensated by a shift of Bint towards smaller values 
for increasing T. Note that both experimental and the- 
oretical uncertainities (precise absolute values of weights 
and fiuctuation corrections, respectively) presently pro- 
hibit a precise quantitative determination of, e.g., the T 
dependence of the parameter t. 
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